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Summary 

A computational method to predict modal reflection coefficients in cylindrical 
ducts has been developed based on the work of Lordi, Homicz, and Rehm, which uses the 
Wiener-Hopf method to account for the boundary conditions at the termination of a thin 
cylindrical pipe. The purpose of this study is to develop a computational routine to 
predict the reflection coefficients of higher order acoustic modes impinging on the 
unflanged termination of a cylindrical duct. 

This effort was conducted under Task Order 5 of the NASA Lewis LET Program, 
Active Noise Control of Aircraft Engines: Feasibility Study , and will be used as part of 
the development of an integrated source noise, acoustic propagation, ANC actuator 
coupling, and control system algorithm simulation. 

The reflection coefficient prediction will be incorporated into an existing 
cylindrical duct modal analysis to account for the reflection of modes from the duct 
termination. This will provide a more accurate, rapid computation design tool for 
evaluating the effect of reflected waves on active noise control systems mounted in the 
duct, as well as providing a tool for the design of acoustic treatment in inlet ducts. As an 
active noise control system design tool, the method can be used preliminary to more 
accurate but more numerically intensive acoustic propagation models such as finite 
element methods. 

The resulting computer program has been shown to give reasonable results, some 
examples of which are presented Reliable data to use for comparison is scarce, so 
complete checkout is difficult, and further checkout is needed over a wider range of 
system parameters. In future efforts the method will be adapted to provide input 
reflection coefficients to the GEAE segmented cylindrical duct modal analysis prediction 
program. 
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Active Control of Fan Noise: Feasibility Study 
Volume 5: Numerical Computation of Acoustic Mode 

Reflection Coefficients for an Unflanged Cylindrical Duct 


1. Introduction 

1.1 Study Assumptions and Objectives 

The purpose of this study was to develop a computational routine to predict the 
reflection coefficients of higher order acoustic modes impinging on the unflanged 
termination of a cylindrical duct. This prediction will be incorporated into the existing 
cylindrical duct modal analysis to account for the reflection of modes from the duct 
termination. This will provide a more accurate, rapid computation design tool for 
evaluating the effect of reflected waves on active noise control (ANC) systems mounted in 
the duct 1 , as well as providing a tool for the design of acoustic treatment in inlet ducts. As 
a treatment design tool, the method can be used preliminary to more accurate but more 
numerically intensive acoustic propagation models such as finite element methods. 

The classical duct modal analysis 2 provides an exact solution to duct propagation 
in treated ducts, accurate to the convergence afforded by the number of modes used in the 
solution expansion. It is limited in the sense of requiring idealized geometry; the duct 
must be cylindrical with constant radius along its axis, and the reflection coefficient is that 
for a cylindrical pipe with infinitely thin walls. Aircraft engine ducts have variable duct 
radius and nacelles with finite thickness and rounded ends. 

The cylindrical duct modal analysis does allow axial discontinuities in acoustic 
treatment properties (segmented treatment), but allows no circumferential nonuniformities 
It calculates one frequency and one spinning mode order at a time, assuming coupling 
among the radial modes but no coupling among spinning modes. 

The modal analysis should be sufficiently accurate to allow evaluation of active 
noise control concepts when adapted as part of a simulation model that includes the 
coupling of the ANC actuators to the duct walls and the ANC feedback control system. If 
required, more accurate prediction and confirmation of ANC or treatment suppression or 
farfield radiation absolute levels cam subsequently be done with a more accurate prediction 
model (such as FEM). 

An objective was to develop a reflection coefficient computational routine that is 
reasonably rapid when coupled to the modal analysis, so that treatment design parametric 
studies involving many cases with different frequencies and spinning mode orders are 
practical and efficient. If only the treatment parameters change, the reflection coefficients 
can be stored and re-used, increasing efficiency for multiple case runs. Twenty years ago, 
when the segmented duct modal analyses were first developed, such a reflection 
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coefficient routine was developed and incorporated into a segmented duct modal analysis 3 , 
but it took so long to compute the coefficients that it was seldom used. 


1.2 The Weiner-Hopf Technique 

The approach in this study is founded upon the prior work of Lordi, Homicz, and 
Rehm 4 , (LHR) particularly as detailed in the second volume of their Air Force contract 
final report 5 . This work is based on the Wiener-Hopf solution for the reflection and 
radiation of waves propagating in a pipe, and it, in turn, is based on the prior work of 
several researchers 6,7,8 . 

The Wiener-Hopf technique 9 is a very complicated mathematical method for 
solving boundary value problems involving semi-infinite waveguides that radiate to the 
farfield. The method transforms the partial differential wave equation and the semi-infinite 
waveguide boundary conditions using spatial Fourier transforms (in wavenumber space) to 
convert the solution to integral equation form. Solution of the integral equation requires a 
special factorization of the kernel, invocation of properties of complex analytic functions, 
and subsequent inversion of the Fourier transform in complex space using contour 
integrals. 

The method, although straightforward in principle, is not unique in the form of 
separating the kernel, and requires application of advanced methods of complex variable 
integral calculus to arrive at a solution that is, more or less, computable. Part of the 
solution manifests itself in closed form, and part of the solution still requires a numerical 
integration. The many steps in the derivation and the invocation of esoteric theorems of 
complex analysis, as well as dealing with Bessel functions and their identities, make the 
entire procedure very error-prone. 

The effort in this development makes no claims of being particularly original. The 
application of the Wiener-Hopf method by Lordi, Homicz, and Rehm to the case of higher 
order spinning modes in a cylindrical pipe was assumed to be derived correctly, and, up to 
the point where they begin to put results in computational form, the derivation was not 
checked step-by-step. LHR present final equations that are almost, but not quite, ready 
for computation. This analysis picks up at a particular stage of the LHR derivation and 
develops a computation that gives reasonable results and runs quickly on an IBM-PC. 
Final debugging and checkout of the computer program is still in progress. Unfortunately, 
reliable data against which to check the computation is scarce. 


1.3 Modal Solution in a Cylindrical Duct 

To define the modal reflection coefficients for wave propagation in a cylindrical 
duct, we present the form of the modal acoustic pressure as given by LHR for the m* 
radial mode of the n* order spinning mode: 
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where 


r = duct radius 
z = axial variable 
t = time 

n = spinning mode order (-00 < n < co) 
m = radial mode number (starting from m = 1 to 00) 

Pnm : (n,m) pressure mode 

= (n,m) pressure mode coefficient 
knm = (n,m) mode axial propagation constant (axial wavenumber) 

J n (x) = n 111 order Bessel function of first kind 
ynm = (n,m) hardwall cylindrical duct eigenvalue 
a = duct radius 

Romj = reflection coefficient for reflection of radial mode m into radial mode j 
co — circular frequency (2nf) 

The complete n** 1 order sp innin g mode solution for the pressure in the duct is the 
superposition of all pressure modes over index m. The cylindrical (inlet) duct geometry is 
shown in Figure 1. 



Figure 1 Duct geometry and coordinate system. 
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Note that the sum over modes involving the reflection coefficient, index j, goes to 
infinity, which implies that all cut-off modes must be included in this summation. A basic 
assumption of this study will be that the reflection coefficient sum can be approximated by 
summing only to the highest cut-on mode order, denoted by mo, so that reflections into 
cut-off modes are neglected (in effect, cut-off mode reflection coefficients are set to zero). 

This will not cause an error in the cut-on mode reflection coefficient calculation, 
and it may be assumed that energy reflected into cut-off modes will attenuate very quickly. 
The reflected cut-off modes may have an effect on the overall segmented duct solution for 
ducts with very short segments at the termination end. Including the computation of cut- 
off mode reflection coefficients should be considered for future program upgrades. 
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2. Basic Equations for the Reflection Coefficient 

Given the spinning mode order, n, the reflection coefficient for the reflection of the 
m th radial mode into the j* radial mode at the termination of an unflanged circular duct, 
based on the LHR analysis, is 



(ka + k^a) 

(ka + k nj a 

)K + (k nm a)K + (k nj a) J n (y^ m ) 

2k nj a | 

k nm a + k„,a) 

^ nm nj j 
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(2-1) 


where 


n = spinning mode order 

m,j = radial mode indices 

Tnm = (n,m) mode eigenvalue 

knm - (n,m) mode axial propagation constant 

a = duct radius 

k = co/c 

c = speed of sound 


The axial propagation constant is given by 


k 


nm a — 



(2-2) 


and the mode eigenvalues yam are determined as the roots of the equation 

Jn(r™)= 0 (2-3) 

These are the eigenvalues for the modes in a hardwall circular duct, determined from the 
boundary condition that requires the normal acoustic velocity (or derivative of pressure) 
to be zero at the wall. 


The K + are defined as functions of argument C& by 

K + (Ca) = ^(ka-Ca)J^(v)H^(v) e^ M(Ca) (2-4) 

for n > 0, where H n (1) is the Hankel function of the first kind, and 


v = v(£a) = y/(ka) 2 ~(£a) 2 (2-5) 

From references on Bessel functions, we can write the Bessel function derivatives and the 
definition of the Hankel function derivative as 
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( 2 - 6 ) 
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V 

Hi (1) (v) = Ji(v) + iY,;(v) 


Yi(v) = -Y„(v)-Y n+1 (v) 

V 


The function M(^a) in Equation (2-4) is given by the integral 
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(2-9) 


where 


k n0 a = ka (2-10) 

otherwise km,, is defined by Equation (2-2). The function M(£a) is a contour integration 
shown by LHR 10 to arise from the Wiener-Hopf kernel factorization procedure. 

The function H(v) is defined by 


where 


Q(v) = arctan 
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( 2 - 11 ) 


Q(0) = 0 

(2-12) 

V( ka ) — ( k nm a ) 

(2-13) 

^( v nm ) ~ ( m - 

(2-14) 


Equation (2-14) defines the branch of the arctangent function to be chosen for each mode, 
and guarantees that £"2(v) is continuous as the arctangent goes through each cycle of n. 
Use of Equation (2-14) necessitates keeping track of the arctangent branch when passing 
from one modal segment to the next in Equation (2-9). 
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The path of the M-integral in Equation (2-9) is shown in Figure 2. The path 
proceeds backwards along the real axis starting at wa = ka and includes all cut-on modes 
before it reaches wa = 0. The index mo denotes the highest order cut-on mode. The 
integral proceeds along the real axis to zero, then heads up the imaginary axis, where it 
picks up all cut-off modes. The parameter Qa is a pole on the real axis (for the reflection 
coefficient calculation Qa = k^a). The integration to high cut-off mode orders should 
eventually converge, since the denominator of the M-integral goes to infinity as i-wa goes 
to infinity. 



Figure 2. M-integral integration path in complex plane (dashed line) 

The parameter Qa takes on different meanings depending upon whether the 
reflection coefficient or the farfield radiation is being calculated. For the reflection 
coefficient calculation, Qa takes on only the values k^a, for all modes m. Thus, it is 
located only at intersections between modal integration segments on the real axis, and the 
poles involve the endpoints of two adjacent integration segments. For the farfield 
radiation calculation, the Qa may occur inside any segment, as well as at the endpoints. 
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3. Behavior of the Integrand of the M-Integral 

The function M(Ca), defined in Equation (2-9), has a direct effect on the value of 
reflection coefficient Rnmj through the K+ factor, defined in Equation (2-4). The 
computation of the integral in M(Ca), which we shall refer to as the M-integral, is the most 
difficult aspect of the entire reflection coefficient calculation. The M-integral can be split 
into two terms: 


2Ca « k T a Q(v) . 

GV = — I f ; ' 2 V ; 2 d(wa) 

71 m=1 (Ca) -(wa) 2 


and 


2Ca 00 k »" a 
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(3-1) 


(3-2) 


where GV is referred to as the variable-Q part of the M-integral and GC as the constant- 
Q part. Both integrals are written as a sum over all cut-on and cut-off modes. 

It is instructive to examine the behavior of the integrand functions of the M- 
integral, particularly the function Q(v), over the real axis from ka down to zero. As a 
sample case, we shall use the following set of parameters: 

Spinning mode order, n = 2 
ka = 12.5 


For these conditions, there are three cut-on modes that have the eigenvalues and 
propagation constants given in Table 1 . : 


Radial mode order m 

Eigenvalue Y 2 m 


1 

3.054237 

12.12112 

2 

6.706133 

10.54883 

3 

9.969468 

7.540437 


Table 1. Test case mode eigenvalues and propagation constants 


We are examining in particular the variable-Q part of the M-integral over the real 
axis for 0 < wa < ka. The integrand for this integral may be written as 

£ Xy) g(v) Q(v) (3 _ 3) 

Ca 2 - wa 2 (Ca - wa)(Ca + wa) Ca - w 
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where Q(v) is the non-singular factor in the integrand, 


Q(v) = 


Q(v) 

Qa + wa 


(3-4) 


We shall examine Q(v), Q(v), and the complete integrand separately, with Qa chosen 
arbitrarily as the second mode wavenumber, Qa = 10.54883. 

The MathCAD program FCM2.MCD, used to carry out the above computations, 
was used to generate the output in Figure 3, and shows what values are obtained when the 
Q-function is calculated directly from the arctangent function definition, using simply the 
arctangent principle value. Notice that as the Bessel function derivatives wind through 
their zero-values, the arctangent has discontinuities when the principle value cycles 
through 7t. The function can be made continuous by adding an appropriate multiple of n, 
depending on how many cycles have been accumulated. 

When Q is re-defined to include these branch cut shifts, the function becomes 
continuous, as illustrated in the FORTRAN-generated plot shown in Figure 4. Notice 
here that the Q-function is continuous at the intersections between mode segments, and 
that the wa-axis is plotted backwards, in the direction of integration. 

Figure 5. shows the non-singular part of the integrand, Q(v), for this case. This 
function is well-behaved within the real axis integration limits, as long as D(v) is properly 
defined. 

Finally, note the behavior of the full integrand in Figure 6 ., with a first-order pole 
at Qa = 10.54883. A special Cauchy principle value integration procedure must be used to 
compute this improper integral containing a singularity. 

The constant-O part of the M-integral differs in that Cl undergoes a discontinuous 
jump in otherwise constant value as the junction of endpoints between two modal 
segments is crossed by the integration. The integrals, however, can be computed in closed 
form within each modal segment, with special attention to the integration across the pole. 

The integrands for the part of the M-integral up the imaginary axis pose no 
additional special problems, and will be considered later. By limiting Qa to real values, we 
are limiting the reflection coefficients to include only cut-on modes. This allows the 
effects of cut-off modes to be included in the computation of each cut-on mode reflection 
coefficient, but omits the reflection coefficients for cut-off modes themselves. This avoids 
the (not too difficult) problem of complex values of Qa, which would lead to poles along 
the imaginary axis. Other than causing programming difficulties with the complex algebra, 
there is no practical reason why cut-off modes could not be included, but this will be left 
for future program upgrades. 
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Figure 5. Function Q(wa' 





4. Computation of the M-Integral 


4.1 Real Axis, Variable Q 

This integral, for each modal segment, can be written 

Otvl 

GMVR(m) = / - Y , d(wa) (4-1) 

(Ca) -(wa) 


where 


k n0 a = ka (4-2) 

and m goes from m = 1 to m = mo, where mo is the highest order cut-on mode. In 
addition, there is the partial segment 


0 fifvl 

GMVR(m 0 +l) = J // d(wa) (4-3) 

k, m „a (Ca) -(wa) 

which completes the full real axis range of 0 < wa < ka. 

Over a “clean” segment that contains no poles at either endpoint, this integral must 
be computed numerically. A Romberg FORTRAN integration subroutine 11 has been 
adopted for this purpose, and it converges rapidly for clean segments. 

The part of the integration that requires special care is the choice of the correct 
branch of fl(v) for the segment. The segment endpoints are obtained from the roots of 
the derivatives of the Bessel function of the first kind (Equation (2-3)). A subroutine to 
compute these roots is provided with the program. 

This Bessel function root extraction program is accurate to about IE-6 in single- 
precision arithmetic, so that when Q(v) is computed close to an endpoint in the numerical 
integration scheme, care must be taken not to slip to the other side of the root, which will 
cause a discontinuity in f2(v) and the integrand. To ensure that a particular numerical 
integration remains entirely within the proper segment, so that Q(v) is computed with the 
correct multiple of n adder, the limits of integration are started and stopped arbitrarily 
0.001 units of wa from the endpoints of the segment. 

This causes a small error in the computation, but avoids a much larger error that 
would result from the Romberg scheme with a discontinuous integrand, not to mention an 
extreme time penalty as the Romberg scheme attempted to iterate through the 
discontinuity. This is a quick, crude, practical solution to the problem, and probably more 
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sophisticated schemes could be developed that would allow closer approach of the integral 
limits to the endpoint roots but still guarantee continuity. 

The integration across the endpoint pole involves two adjacent segments that must 
be considered together. A pre-integration sorting scheme identifies each segment as one 
of the following types: 


Type 1 - Clean segment, no poles 
Type 2 - Pole internal to endpoints 

Type 3 - Pole at upper integration limit endpoint (lower wa- value) 

Type 4 - Pole at lower integration limit endpoint (higher wa-value) 

The integral across the pole is considered in the Cauchy principle value sense 12 , which 
makes the basic assumption that the integral exists. We know that the Q(v) factor in the 
integrand is continuous, so that the only problem arises from the l/(£a - wa) singularity 

As an illustration of the application of the Cauchy principle value theorem to the 
integration across a pole, consider the simpler form 


I 



dx 

x - 2 


(4-4) 


This can be computed as 

= g * 0 {M* - 2 l]?" e + N* - 2 fl2+e } = 
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1 = L 
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dx 
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(4-5) 


The terms in e can be cancelled before taking the limit, giving, simply 

1= f-^- = ln3 = 1.10 (4-6) 

J x - 2 

In the M-integral, the function Q(v) in the integrand prevents integration across the 
pole in closed form, so that an approximate method must be used. Consider the 
integration across the pole between arbitrary limits ±s on either side of the pole, as shown 
in Figure 7. The integral is written as 

^ a— e Ofvl 

SG = J -^-d(wa) (4-7) 

Ca +e Ca-wa 
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Figure 7. Integration across singularity on real axis 

If s is chosen sufficiently small, we can approximate Q(v) by a linear function, as 
illustrated in Figure 8 and written as 

Q(x) = Djwa + Dq (4-8) 



Figure 8. Linear fit to Q-function around pole 
Then we can write 


-V DlWa+D ° d(wa) = 

Ca+s wa "£a 


^ a f G Djwa t- 

f — l - d(wa) + 

Ca+e wa - C a l 
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When carried out in the principle value sense, this can be shown to give simply 


SG = 2D ie (4-10) 

In the computer program, this contribution from the integral around the pole is 
added to the numerical integration for segment j taken up to C , a + e, and the numerical 
integration for segment j+1 includes only the contribution from the integral starting at 
Cfl. - 8 . 


The coefficients D 0 and Di are obtained from a simple 3 -point linear regression 
curve fit to the values of Q(v) at the points 

P 1 : wai = Qa - s 

P2: wa 2 = Qa 

P3 : wa 3 = Qa e 

as shown in Figure 8. An example of the accuracy of the linear assumption using s = 0.05 
is shown in Figure 9. 

Examination of a few cases has indicated that s = 0.05 is a reasonable assumption, 
introducing only small errors in the overall integral. The trade-off is in computation time, 
since the Romberg numerical integration scheme requires many more iterations to 
converge to the specified tolerance the closer the endpoint comes to (i.e., starts going up) 
the pole. A diagnostic program is available to examine the details of the variable-fi 
integration along the real axis. 
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Figure 9. Linear fit of function Q(wa) around pole for m = 2, C,a = 10.54883, ka = 
12.5, with e = 0.05. 


4.2 Real Axis, Constant H 

The integral along the real axis with the constant-Q numerator in the integrand can 
be written 


GMCR(m) = J 


Q ( v n,m-l) 
(C a ) 2 - (wa) 2 


d(wa) 


(4-11) 


where Equation (4-2) still applies and there is an integral analogous to Equation (4-3) over 
the partial segment from k„, m0 a to zero. The integral for GMCR can be computed in 
closed form, with special attention to integrating across any singularities. 

Part of this integral, it will turn out, is required in a special form to be used to 
resolve an indeterminate in the reflection coefficient K + function. The form needed for this 
integral is solved by LHR (Reference 4, p. 48, Eqn (A-37)), and is given as 


GMCR(m) = X In 
m=l 


k nm a + C a 
k nm a - C a 


+ i g(v n ) + n(v n m )j 


(4-12) 
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where mi is the mode whose axial wavenumber is just greater than or equal to Qa. This 
integral includes the 2£a/7t multiplication factor that is not yet included in the other 
integrals. 

To derive the integral in this form, LHR multiply the integrand in Equation (4-11) 
by 2C,a and expand as the sum 


2Ca = 1 | 1 

i^a 2 - wa 2 C a + wa Qa - wa 


(4-13) 


They then operate on each of the resulting two integrals separately, solving each in closed 
form to give GMCR. 

The summation of natural log term, when substituted into the exponential in 
M(Ca), (see Equation (2-4)), becomes the factor 


exp 


m o 

I 

In 

^nm a + C a 

m o 

= n 

knm a + C a 

vm=l 


_ k n m a - C a _, 

' m=l 

knma-^a 


(4-14) 


One of the factors in this multiple product will be combined with the factor in the 

reflection coefficient to resolve an indeterminate form. The remaining two constant terms 
are purely imaginary, and will affect only the phase of the reflection coefficient when 
substituted into the M-integral exponential. 


4.3 Imaginary Axis, Variable Q 

For the integration over modes up the imaginary axis, we make the transformation 
of variable 

wa = iy (4-15) 

so that y is a real variable of integration. Inserting this into the integral for each modal 
segment gives 

GMVI0) = iT -^Idy (4-16) 

aC a 2 +y 2 

where index j implies a cutoff mode order greater than mo and 

i K nj a = k n j a = i-^Ynj "(ka) 2 (4-17) 
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making K n j real, since, for cutoff modes, Ynj > ka. 

The first value of j in Equation (4-16) will be j = mo+2, that is, the second mode 
order above cutoff. We must also include the partial segment integral starting from y = 0, 

^n.iriO+A OfvfvYl 

GMVT(m 0 + 1) = i j ; ^ dy (4-18) 

0 & + T 

In the sum of integrals over all real and imaginary mode segments, GMVR(mo+l) from 
Equation (4-3) combines with GMVI(mo+l) to give the total contribution of mode mo+1 
to the M-integral. 

When transformed to the y-variable, the argument of the Bessel function in fi(v) 
becomes 


v(wa) = v(y) = ^/(ka) 2 - (iy) 2 = -y/(ka) 2 + y 2 (4-19) 

which is real. Thus, the entire integrand of Equation (4-16) is real, and GMVI(j) is purely 
imaginar y This indicates that integral up the imaginary axis contributes only to the phase 
of the reflection coefficient, through the M-integral exponential. 

It is assumed without proof that the sum of integrals up the imaginary axis 
converges to a finite value, and that only a finite number of terms need be included to 
obtain sufficient accuracy. This assumption will be examined later in test case 
computations. 

Since D(v) is variable, the integration in Equations (4-16) and (4-18) must be 
computed numerically. This is accomplished using the Romberg method, which can be 
carried out as a real integral. All integrals are of the “clean segment” type, with no 
singularities, since Qa is real. 

If it were desired to compute the reflection coefficient for a cut-off mode, such 
that C,a became imaginar y, a method similar to the method of integration through a 
singularity on the real axis would have to be used on the imaginary axis. To reiterate, 
although reflection coefficients for scattering into or among cut-off modes are not 
calculated, the effect of the scattering into cutoff modes on the reflection coefficient of the 
cut-on modes (a phase shift effect) is included. 


4.4 Imaginary Axis, Constant Q 

This integral, for each mode segment, can be written 
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K -) a n(v„i_,) 

GMCI(j) = i ( -V^-fdy (4-20) 

c* *r 

where the transformation of Equation (4-15) has been used to make the integral itself real. 
As in the case of GMVI, the integral over the partial segment from y = 0 to y = Kn, m o + ia 
must also be included. 

There are no singularities in this integral, and, since the numerator is a constant 
over the mode segment, it can be evaluated in closed form as 

M v ni-l)r ( K ni a> l K n_ i — l a> l 

GMCI(j) = i — - — - — - arctan — — I - arctan — — (4-21) 

C a \_ v Ca ) Ca J 

LHR integrate this part using the same expansion of the denominator into two 
terms that they used to find GMCR, leading to a log sum solution. The form in Equation 
(4-21) is equivalent to the LHR form, and should be slightly easier to compute. No 
problem should be encountered by taking the principle value of the arctan function in 
Equation (4-21). 


4.5 Final Form of the M-Integral 

After shifting the multiple product factor in GMCR from the exponential of the M- 
integral to the multiplier of the exponential, the remainder of the M-integral can be written 
as 


2Ca m o 

M(Ca) = — \ X GMVR(m) + GMVR(m 0 +1) 
71 [m=l 

2Ca f Max ] 

+ — \ X GMVI(j) + GMVI(mo +1) 

n [j=m 0 +2 J 

2Ca f Max 1 

- — \ X GMCI(j) + GMCI(m 0 + 1) 


i [ i X V n,m 0 ) + n,m 1 )] 


% [j= m o + 2 


where expressions for the integrals are derived above. Integrals that contain a “V” in the 
name must be evaluated numerically, while integrals with names containing a “C” are 
given by closed form expressions. The required procedures for integrating across 
singularities must be observed for computing the GMVR integrals. The sums of integrals 
up the imaginary axis are continued to some maximum value Max, and this must be 
suitably high to provide an acceptable convergence tolerance. 
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5. Computation of Reflection Coefficient 
5.1 Final Computational Form 

When the multiple product derived from the M-integral integration of the constant- 
Q part of the integration along the real axis is combined with the K+(£a) function from 
Equation (2-4), we obtain a new form: 


K + fta) = ^n(ka-Ca)Ji(v)Hi (1) (v) 


lll 0 

n 

m=l _ 


k nm a +Ca 
k nm a -Ca 


e iM(Ca) 


(5-1) 


where now M(Ca) is the revised form given by Equation (4-22). Combining the radicals 
and rearranging slightly gives the square-root factor 


:<ka -Ca) J ° (V) H^(v)(k nj +Ca) U 
knj a 0> a m=l 

m*j 

where the j* factor in the multiple product has been extracted and its denominator 
associated with J^(v) . 



k nm a + C a 

L k nm a — C a 


(5-2) 


As Ca => k nj a, this ratio becomes an indeterminate form (0/0), since v n j = y„j is a 
root of the Bessel function derivative. Using L’Hospital’s procedure for indeterminate 
forms (taking the derivative of the numerator and the denominator with respect to Ca 
before taking the limit), we can show that 


l 

;a-»k n ja 


J n ( v (C a )) 
k nj a -£a 


= i|^(n-Ynj)jn(Ynj) 
Ynj 


(5-3) 


Since Jn(Ynl) = 0 , we have 


H^Ynj^iY^Ynj) 


Using the Bessel function identity 


J n (x)YA(x)-JK(x)Y n (x) = — 

7DC 


(5-4) 


(5-5) 


we have 
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2i 


(5-6) 


* Jn(7 nj)^n(Y nj) 


^nj 


Using Equations (5-2) through (5-6), the product of the K + functions can finally be written 
as 


K.(k m ,a)K.(k„ 11 a) = 


4 ‘ k r a T^~ i/(ka - k^aXka - k^aXn 2 - T» X° ; - t! )n„n 

y nfi y nm 


(5-7) 


where 


nip 


n„ = n 


j=l 

j*m 


k DJ a+Ca 

k nj a-C a 


(5-8) 


and 


flip 

n,=n 


j=l 


M-Ca 


(5-9) 


This gives the final computational form for the reflection coefficient for the m 11 ' 
radial mode reflecting into the p* radial mode: 


o 2knin Jn(y ran^) 

nnv y 2 (k a + k a)j fy a) 

I nm\ ran n|i“/ n \‘ n»» / 


(5-10) 


((ka) J -(k^a) ! )((ka) ! -(k m a) 2 )^Xf-n„n„ 


M(k^a) 


where M(k n ^a) is obtained from Equation (4-22). 


5.2 Sample Calculations 

An edited computer printout follows of the computation of the reflection 
coefficient for a test case, denoted Case #1, with n = 2, ka = 12.5, and s = 0.05, where e is 
the distance around the singularity used for the linear fit. Diagnostic results for the 
intermediate computations of the integrals are included in the Case #1 printout. 
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Note that for this value of ka, no modes are near cut-off, so that all reflection 
coefficient magnitudes are small, the largest being the R 13 cross-term with a magnitude of 
0.132. Note that from the results of the imaginary axis integrations that the last term for 
the 10 th cut-off mode is about 3-5% of the sum for the constant-^ or variable-Q parts. 


PROGRAM CDRC1 TEST CASE #1 

CALCULATES MODAL REFLECTION COEFFICIENTS 
USING HOMICZ & LORDI M-INTEGRAL FORM 

SPINNING MODE ORDER M = 2 

ka = 12.50000000 

NUMBER OF CUT-ON MODES = 3 

NUMBER OF MODES IN CALC = 13 

EPSILON ABOUT POLE - .050000 

EIGENVALUES AND PROPAGATION CONSTANTS 
ROOTA ARE INTEGRATION SEGMENT ENDPOINTS 


MODE 

EV 

PKA REAL 

PKA IMAG 

ROOTA 

1 

3.05423700 

12.12112000 

.00000000 

12.12112000 

2 

6.70613300 

10.54883000 

.00000000 

10.54883000 

3 

9.96946800 

7.54053700 

.00000000 

7.54053700 

4 

13.17037000 

-.00000018 

4.14833500 

.00000000 

5 

16.34752000 

-.00000046 

10.53525000 

4.14833400 

6 

19.51291000 

-.00000065 

14.98345000 

10.53525000 

7 

22 . 67158000 

-.00000083 

18.91430000 

14.98345000 

8 

25.82604000 

-.00000099 

2 2.59943000 

18.91430000 

9 

28.97767000 

-.00000114 

26.14298000 

22.59943000 

10 

32.12733000 

-.00000129 

29.59586000 

26.14298000 

11 

35.27554000 

-.00000144 

32.98656000 

29.59586000 

12 

38.42265000 

-.00000159 

36.33250000 

32.98656000 

13 

41.56894000 

-.00000173 

39.64501000 

36.33250000 

14 

.00000000 

.00000000 

.00000000 

39.64500000 


REAL AXIS INTEGRATIONS: 

GV2 ARE SEGMENT INTEGRALS FOR VARIABLE OMEGA PART 
JS = SEGMENT #: ISX = OMEGA MULTIPLIER: ITP = SEGMENT TYPE 

FOMM0 AND FOMMZ ARE TERMS IN CONSTANT OMEGA INTEGRATION 
JM = RADIAL MODE NUMBER: GMR * TOTAL REAL AXIS INTEGRAL 

JS ISX ITP GV2 

103 -3. 199366 0E-02 

2 14 -1 . 4 871500E-01 

321 -2 . 3233030E-01 

431 -4.7253740E-01 

JM= 1 FOMM0= 6 . 2832600E+00 FOMMZ= 2 . 5844570E-04 

REAL AXIS INTEGRAL, JM= 1 GMR= - 6 . 8335 92 0E+00+ 6 . 2 835 190E+00*i 

JS ISX ITP GV2 

101 -2.2201310E-03 

2 1 3 3 . 7 623970E-01 

3 2 4 -8 . 1834 660E-01 

431 -6. 6436000E-01 

JM= 2 FOMM0= 6 . 2832600E+00 FOMMZ = 3 . 1 4 174 10E+00 
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REAL AXIS INTEGRAL, JM= 2 GMR= -7 . 4454 900E+00+ 9. 4250010E+00*i 

JS ISX ITP GV2 

1 0 1 -9.2 672090E-04 

211 4 . 1502940E-02 

323 1. 4690930E+00 

434 -2 . 7157840E+00 

JM- 3 FOMMO= 6 . 2 832 600E+00 FOMMZ= 6 . 2 832 600E+00 

REAL AXIS INTEGRAL , JM= 3 GMR- -5 . 78 99020E+00+ 1 . 2566520E+01*i 

IMAGINARY AXIS INTEGRALS FOR 10 CUTOFF MODES 
XI & X2 ARE SEGMENT ENDPOINTS ON IMAGINARY AXIS 
GCI IS CONSTANT OMEGA PART; GVI = VARIABLE OMEGA PART 
NOTE THAT GCI AND GVI ARE IMAGINARY NUMBERS 


JS 

ISX 

XI 

X2 

GCI 

GVI 

NST 

1 

3 

.0000000 

4.1483340 

2 . 38 42440E-01 

2 . 4418080E-01 

2 

2 

4 

4.1483340 

10.5352500 

2.9995 660E-01 

3 . 4049100E-01 

2 

3 

5 

10.5352500 

14.9834500 

1 . 8 152870E-01 

2 . 0208300E-01 

1 

4 

6 

14.9834500 

18.9143000 

1 . 4290520E-01 

1 . 5 618500E-0 1 

1 

5 

7 

18.9143000 

22.5994300 

1.2069470E-01 

1 . 3015 610E-01 

1 

6 

8 

22.5994300 

26. 1429800 

1.0551070E-01 

1 . 12 65 110E-01 

1 

7 

9 

26.1429800 

29.5958600 

9 . 4195800E-02 

9 . 9799390E-02 

1 

8 

10 

29.5958600 

32.9865600 

8 . 5317340E-02 

8 . 9842860E-02 

1 

9 

11 

32.9865600 

36.3325000 

7. 8105830E-02 

8 . 1842490E-02 

1 

10 

12 

36.3325000 

39.6450000 

7 . 2101210E-02 

7 . 5241620E-02 

1 


INTEGRAL TOTALS ( JM, GVI , GCI , GMISUM) 

1 1 . 5324740E+00 1 . 4 1874 10E+00 2 . 2773180E+01 


M- INTEGRAL FOR MODE JM= 1 

GM= -6 . 8335920E+00 + 2 . 905 6700E+01* i 


JS 

ISX 

XI 

X2 

GCI 

GVI 

NST 

1 

3 

.0000000 

4.1483340 

3 . 1129010E-01 

3 . 1873770E-01 

2 

2 

4 

4.1483340 

10.5352500 

3 . 6638410E-01 

4 . 1483920E-01 

1 

3 

5 

10.5352500 

14.9834500 

2 . 05 63030E-01 

2 . 2870910E-01 

1 

4 

6 

14.9834500 

18.9143000 

1 . 55 87390E-01 

1 . 7028760E-01 

1 

5 

7 

18.9143000 

22.5994300 

1 . 2871 610E-01 

1 . 3877520E-01 

1 

6 

8 

22.5994300 

26.1429800 

1.1088770E-01 

1 . 1837650E-01 

1 

7 

9 

26.1429800 

29.5958600 

9 . 8001730E-02 

1.0382320E-01 

1 

8 

10 

29.5958600 

32.9865600 

8.8121140E-02 

9 . 2790350E-02 

1 

9 

11 

32.9865600 

36.3325000 

8 . 0236280E-02 

8 . 4071740E-02 

1 

10 

12 

36.3325000 

39.6450000 

7.3760720E-02 

7.6971370E-02 

1 


INTEGRAL TOTALS {JM, GVI , GCI , GMISUM) 

2 1 . 7473820E+00 1 . 618 902 0E+00 2 . 2 6065 90E+01 


M- INTEGRAL FOR MODE JM= 2 

GM= -7 . 4454 900E+00 + 3 . 20315 90E+01 *i 


JS 

ISX 

XI 

X2 

GCI 

GVI 

NST 

1 

3 

.0000000 

4.1483340 

5 . 8457050E-01 

5 . 9809270E-01 

2 

2 

4 

4.1483340 

10.5352500 

5.5824050E-01 

6 . 2806500E-01 

2 

3 

5 

10.5352500 

14.9834500 

2 . 5824970E-01 

2 . 8668140E-01 

2 

4 

6 

14.9834500 

18.9143000 

1 . 8100920E-01 

1 . 9758640E-01 

2 

5 

7 

18.9143000 

22.5994300 

1. 4326990E-01 

1 . 5440400E-01 

1 

6 

8 

22.5994300 

26. 1429800 

1 . 2024990E-01 

1.2834160E-01 

1 

7 

9 

26.1429800 

29.5958600 

1 . 0444730E-01 

1. 1063630E-01 

1 

8 

10 

29.5958600 

32.986560 0 

9 . 2777270E-02 

9 . 7684410E-02 

1 

9 

11 

32.9865600 

36.3325000 

8 . 3723440E-02 

8 . 7720260E-02 

1 

10 

12 

36.3325000 

39.6450000 

7 . 64 47380E-02 

7 . 9771590E-02 

1 


INTEGRAL TOTALS (JM, GVI , GCI , GMISUM) 

3 2 . 3689840E+00 2 . 20298 60E+00 2 . 1947530E+01 
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M- INTEGRAL FOR MODE JM= 3 

GM= -5.7899020E+00 + 3 . 45 14050E+01*i 


REFLECTION COEFFICIENTS - PHASE IN DEGREES 


JM 

JL 

REAL-RC 

IMAG-RC 

RCMAG 

RCPHS 

1 

1 

4 . 7273 92 0E-02 

4 . 698 8 270E-02 

6 . 6653740E-02 

44.8264 

1 

2 

5.842938 0E-02 

4 . 911 6320E-02 

7 . 63308 90E-02 

40.0507 

1 

3 

-1.2287330E-01 

- 4 . 75 8 0420E-02 

1 . 3176390E-01 

-158.8320 

2 

1 

3.3685 690E-02 

2 . 831 6530E-02 

4 . 4 006270E-02 

40.0507 

2 

2 

-4.1341 0 60E-02 

-2.9244 140E-02 

5 . 0638940E-02 

-144.7249 

2 

3 

-8 . 4828000E-02 

-2 . 4 953 930E-02 

8 . 8422210E-02 

-163.6076 

3 

1 

-3 . 52 155 00E-02 

-1.3 63 65 60E-02 

3 . 7763570E-02 

-158.8320 

3 

2 

-4.2 1 69790E-02 

- 1 . 2 405 130E-02 

4 . 395 655 0E-02 

-163.6076 

3 

3 

7 . 88 62390E-02 

-3 . 4299810E-03 

7 . 8936950E-02 

-2.4904 


Case #2 below is the same as Case #1 except that s is reduced to 0.02. Comparing 
the reflection coefficients for this case with Case #1 indicates very small differences in the 
reflection coefficient magnitudes and no difference in the phases. This indicates that, at 
least for n = 2, s = 0.05 is a reasonably good assumption, and the total reflection 
coefficient computation at £ = 0.05 takes only a few seconds on a 33MHz IBM-PC 486. 


PROGRAM CDRC1 TEST CASE 2 

CALCULATES MODAL REFLECTION COEFFICIENTS 
USING HOMICZ 4 LORD I M- INTEGRAL FORM 

SPINNING MODE ORDER M = 2 

ka = 12.50000000 

NUMBER OF CUT-ON MODES = 3 

NUMBER OF MODES IN CALC = 13 

EPSILON ABOUT POLE = .020000 


EIGENVALUES AND PROPAGATION CONSTANTS 


MODE 

EV 

PKARL 

PKAIM 

ROOTA 

1 

3.05423700 

12.12112000 

.00000000 

12.12112000 

2 

6.70613300 

10.54883000 

.00000000 

10.54883000 

3 

9.96946800 

7.54053700 

.00000000 

7.54053700 

4 

13.17037000 

-.00000018 

4.14833500 

.00000000 

5 

16.34752000 

-.00000046 

10.53525000 

4.14833400 

6 

19.51291000 

-.00000065 

14.98345000 

10.53525000 

7 

22.67158000 

-.00000083 

18.91430000 

14.98345000 

8 

25.82604000 

-.00000099 

22.59943000 

18.91430000 

9 

28.97767000 

-.00000114 

26.14298000 

22.59943000 

10 

32.12733000 

-.00000129 

29.59586000 

26.14298000 

11 

35.27554000 

-.00000144 

32.98656000 

29.59586000 

12 

38 . 42265000 

-.00000159 

36.33250000 

32.98656000 

13 

41.56894000 

-.00000173 

39.64501000 

36.33250000 

14 

.00000000 

.00000000 

.00000000 

39.64500000 
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REFLECTION COEFFICIENTS 


JM 

JL 

RL-RC 

IM-RC 

RCMAG 

RCPHS 

1 

1 

4 . 728 6910E-02 

4 . 700 1180E-02 

6.6672050E-02 

44.8264 

1 

2 

5 . 845Q030E-02 

4 . 9133 670E-02 

7 . 6357860E-02 

40. 0507 

1 

3 

-1 . 2288 45 0E-01 

-4 . 7584760E-02 

1 . 3177600E-01 

-158. 8320 

2 

1 

3 . 3 6975 90E-02 

2 . 832 6530E-02 

4 . 4021820E-02 

40.0507 

2 

2 

-4 . 135 8 900E-02 

-2.925 6770E-02 

5 . 066Q800E-02 

-144.7249 

2 

3 

-8.4 8 42 4 00E-02 

-2.4958 170E-02 

8 . 8437220E-02 

-163. 6076 

3 

1 

-3 . 5218710E-02 

-1 . 3637800E-02 

3 . 7767010E-02 

-158.8320 

3 

2 

-4 . 2 176950E-02 

-1 . 24 07230E-02 

4 . 3964010E-02 

-163.6076 

3 

3 

7 . 8855 130E-02 

-3 . 42 9 6650E-03 

7 . 8929690E-02 

-2. 4904 


The reflection coefficient Rjm,xl is the factor by which the JM* forward-travelling 
radial mode is multiplied to give the contribution of the JM forward mode to the 
backward-travelling JL mode. The magnitude of the complex value Rjm,jl gives the factor 
by which the amplitude of the component of the reflected mode decreases, and the phase 
of Rjm,jl is gives the phase shift of the reflected mode component. 

Figure 10 shows the computation of the R 33 reflection coefficient magnitude and 
phase as a function of ka for frequencies just above cut-off of the n = 2, j = 3 mode to a ka 
value just below cut-on for the j = 4 mode. The reflection coefficient appears to be 
tending toward a value of 1.0 in magnitude and -180° in phase at the cut-off ka- value. As 
ka approaches the cut-off ka, the value of e must be decreased to keep it between ka and 
the cut-off point. 
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6. Conclusions and Recommendations 


A computational method to predict modal reflection coefficients in cylindrical 
ducts has been developed and gives reasonable results. Reliable data to use for 
comparison is scarce, so complete checkout is difficult. Further checkout is needed over a 
wider range of system parameters. One suggestion is to compare the computations versus 
a finite element prediction model. 

The method will be adapted as a subroutine to the GEAE segmented cylindrical 
duct modal analysis. 

A near-term extension would be to apply the analysis to the prediction of the 
complex farfield directivity pattern, which can be done directly. The M-integral in the 
farfield case should be the same, only the algebraic factor is different. There may be other 
indeterminate forms to be resolved. 

Upgrades that should be considered are to include the effect of mean flow in the 
duct and to examine the difficulty of extending the analysis to annular ducts. 

The computation could be improved by adjusting the starting value of the 
integration along the real axis based on the behavior of the fi-function for wa near ka. In 
this case, Y’„ goes to infinity but Cl goes to zero. Currently, the arbitrary value of 
wa = ka - 0.05 is used, but a better method to determine the starting value has been 
developed that could be incorporated into the numerical integration subroutine. 

The analysis could be upgraded to double precision, particularly the determination 
of the Bessel function roots, but there is probably more error in the integral 
approximations and the assumption of an idealized unflanged pipe than this would 
warrant. Numerical precisions in the computation are currently fairly well “balanced”. 

The method could be extended to include a variable number of cut-off modes in 
the integral up the imaginary axis, so that convergence properties could be examined. The 
computation could also easily be extended to predict the reflection coefficient for cut-off 
modes. 


28 



7. References 


Kraft, R. E., and Kontos, K. B., “Theoretical Implications of Active Noise Control 
for Turbofan Engines”, AIAA 93-4355, October, 1993. 

Rice, E. J., “Propagation of Waves in an Acoustically Lined Duct with a Mean 
Flow”, NASA SP-207, July, 1969, pp. 345-355. 

Zorumski, W. E., “Generalized Radiation Impedances and Reflection Coefficients 
of Circular and Annular Ducts”, J. Acoustical Soc. Amer, V54, No. 6, 1973, pp. 
1667-1673. 

Lordi, J. A., Homicz, G. F , and Rehm, R. G., “Effects of Finite Duct Length and 
Blade Chord on Noise Generation by a Rotating Blade Row”, AIAA 74-555, 
April, 1974. 

Lordi, J. A., Homicz, G. F., and Rehm, R. G., “Theoretical Studies on Fan Noise 
Generation by a Transonic Compressor Blade Row”, Calspan Corporation, Report 
No. AI-2836-A-2, August, 1973. 

Noble, B., Methods Based on the Wiener-Hopf Technique , Pergamon Press, 1 958. 
Levine, H. and Schwinger, J., “On the Radiation of Sound from an Unflanged 
Circular Pipe”, Phys. Rev., 73, 1948, pp. 383+. 

Weinstein, L. A., Theory of Diffraction and the Factorization Method 
(Generalized Wiener-Hopf Technique), Golem Press, 1969. 

Crighton, D. G., Dowling, A. P., Ffowcs Williams, J. E., Heckl, M., and 
Leppington, F. G., Modem Methods in Analytical Acoustics: Lecture Notes, 
Springer- Verlag, 1992 (Chapter 5.). 

Lordi, J. A., Homicz, G. F., and Rehm, R. G., “Theoretical Studies on Fan Noise 
Generation by a Transonic Compressor Blade Row”, Calspan Corporation, Report 
No. AI-2836-A-2, August, 1973, p. 44. 

Press, W. H., Teukolsky, S. A., Vetterling, W. T., and Flannery, B. P., Numerical 
Recipes in FORTRAN: The Art of Scientific Computing, Cambridge University 
Press, 1992, p. 134ff. 

Jeffrey, A., Complex Analysis and Applications, CRC Press, 1992, p. 382. 


29 



REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No. 0704-0168 


Public reporting burden for this collection of information is estimated to average 1 hour per response. Including the time for reviewing Instructions, searching existing data sources 
gatheringand maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this 
collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson 
Davis Highway, Suite 1204, Arlington. VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704-0188), Washington. DC 20503. 


1. AGENCY USE ONLY (Leave blank) 


2. REPORT DATE 

September 1996 


4. TITLE AND SUBTITLE 


Active Control of Fan Noise: Feasibility Study 

Volume 5: Numerical Computation of Acoustic Mode Reflection Coefficients 
for an Unflanged Cylindrical Duct 



3. REPORT TYPE AND DATES COVERED 

Final Contractor Report 


5. FUNDING NUMBERS 


WU-538-03-1 1 
C-NAS3-26617 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESSES) 

General Electric Aircraft Engines 
1 Neumann Way 
RO. Box 156301 
Cincinnati, Ohio 45214-6301 


9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135-3191 


8. PERFORMING ORGANIZATION 
REPORT NUMBER 


E-10379 


10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 


NASA CR-198513 


11. SUPPLEMENTARY NOTES 


Project Manager, Laurence J. Heidelberg, Propulsion Systems Division, NASA Lewis Research Center, organization code 
2770,(216)433-3859. 


12a. DISTRIBUTION/ AVAILABILITY STATEMENT 


12b. DISTRIBUTION CODE 


Unclassified -Unlimited 
Subject Categories 07 and 71 

This publication is available from the NASA Center for AeroSpace Information, (301) 621-0390. 


13. ABSTRACT (Maximum 200 words) 

A computational method to predict modal reflection coefficients in cylindrical ducts has been developed based on the work 
of Homicz, Lordi, and Rehm, which uses the Wiener-Hopf method to account for the boundary conditions at the termination 
of a thin cylindrical pipe. The purpose of this study is to develop a computational routine to predict the reflection coefficients 
of higher order acoustic modes impinging on the unflanged termination of a cylindrical duct. This effort was conducted under 
Task Order 5 of the NASA Lewis LET Program, Active Noise Control of Aircraft Engines: Feasibility Study , and will be used 
as part of the development of an integrated source noise, acoustic propagation, ANC actuator coupling, and control system 
algorithm simulation. The reflection coefficient prediction will be incorporated into an existing cylindrical duct modal analy- 
sis to account for the reflection of modes from the duct termination. This will provide a more accurate, rapid computation 
design tool for evaluating the effect of reflected waves on active noise control systems mounted in the duct, as well as provid- 
ing a tool for the design of acoustic treatment in inlet ducts. As an active noise control system design tool, the method can be 
used preliminary to more accurate but more numerically intensive acoustic propagation models such as finite element meth- 
ods. The resulting computer program has been shown to give reasonable results, some examples of which are presented. 
Reliable data to use for comparison is scarce, so complete checkout is difficult, and further checkout is needed over a wider 
range of system parameters. In future efforts the method will be adapted as a subroutine to the GEAE segmented cylindrical 
duct modal analysis program. 


14. SUBJECT TERMS 

Acoustic ducts; Fan noise; Ducted fans 


17. SECURITY CLASSIFICATION 18. SECURITY CLASSIFICATION 19. SECURITY CLASSIFICATION 
OF REPORT OF THIS PAGE OF ABSTRACT 

Unclassified Unclassified Unclassified 


15. NUMBER OF PAGES 

31 


16. PRICE CODE 

Ao3 


20. LIMITATION OF ABSTRACT 


Standard Form 298 (Rev. 2-89) 
Prescribed by ANSI Std. Z39-18 
298-102 
















